Many-Body Interactions of Neutrinos with Nuclei - Observables 
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Background: The total inclusive cross sections obtained for quasielastic (QE) scattering in the Mini Booster 

Neutrino Experiment (MiniBooNE) are significantly larger than those calculated by all models based on the 

impulse approximation and using the world average value for the axial mass of Ma ~ 1 GeV. This discrepancy 

has led to various, quite different explanations in terms of increased axial masses, changes in the functional 

form of the axial form factor, increased vector strength in nuclei, and initial two-particle interactions. This is 

disconcerting since the neutrino energy reconstruction depends on the reaction mechanism. 

Purpose: We investigate whether exclusive observables, such as nucleon knock-out, can be used to distinguish 

between the various proposed reaction mechanisms. We determine the influence of 2p-2h excitations on the 

neutrino energy reconstruction. 

Method: We extend the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) model by explicitly incorporating 

initial 2p-2h excitations. 

Results: We calculate inclusive cross sections and numbers and spectra of knock-out nucleons and show their 

sensitivity to the presence of 2p-2h initial excitations. We also discuss the influence of 2p-2h excitations on the 

neutrino energy reconstruction. 

Conclusions: Inclusive double-differential cross sections, depending only on muon variables, are fairly insensitive 

to the reaction mechanism. 2p-2h excitations lead to an increase in the number n of knock-out nucleons for n > 2, 

while only n = 1 knock-out remains a clean signal of true QE scattering. The spectra of knock-out nucleons are 

also changed, but their shape is hardly affected. In the energy reconstruction, 2p-2h interactions as well as A 

excitations lead to a downward shift of the reconstructed energy; this effect 2p-2h excitations disappears at higher 

energies because the 2p-2h influence is spread out over a wider energy range. 

PACS numbers: 25.30.Pt,13.15.+g,14.60.Pq 



I. INTRODUCTION 

One of the essential ingredients for any extraction of 
the neutrino masses from oscillation experiments is the 
neutrino energy. This energy is not known a priori in 
present-day experiments, because neutrino beams are 
quite broad in energy due to their production mecha- 
nisms. While at higher energies calorimetric methods 
may play a role, at lower energies (a few hundred MeV 
to a few GeV) quasielastic (QE) scattering has been used 
to determine the incoming neutrino energy on an event - 
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by-event basis. This method relies on an identification of 
the reaction mechanism (interaction of the neutrino with 
a single nucleon). It also relies on the use of quasifrcc 
kinematics that describes neutrino scattering on a single, 
free nucleon at rest, thus neglecting any Fermi-motion 
effects; binding is taken into account only by a constant 
removal energy. 

In theoretical calculations the QE cross section is de- 
termined by an interplay of vector and axial couplings 
with their corresponding form factors. The vector cou- 
plings can be rather well determined from electron scat- 
tering experiments on the nucleon that work at a fixed 
energy and permit one to determine the relevant kine- 



matic parameters, i.e. energy and momentum transfer, 
in each event. The corresponding form factors have been 
shown to have a complicated, non-dipole form pQ. For 
the axial couplings the situation is less well determined. 
Here the data come from electro-pion production and 
older neutrino data on the nucleon or deuterium with 
large uncertainties. They have been analyzed by mak- 
ing a dipole ansatz and then extracting the axial mass 
from a fit to data. The world average for the axial mass 
parameter is found to be Ma — 1.026 GeV [2]. 

It came, therefore, as a surprise when the Mini Booster 
Neutrino Experiment (MiniBooNE) at Fermilab pub- 
lished its results on QE scattering. The analyses of both 
charged current (CC) and neutral current (NC) high- 
statistics QE events showed a clear excess of cross section 
over that expected for QE scattering from a Fermi-gas 
model [SHS]. A similar result had been obtained by the 
K2K experiment [6] that worked with a neutrino flux 
peaked at a slightly higher neutrino energy (« 1 GeV) 
than the MiniBooNE experiment (« 0.7 GeV). In both 
cases the flux distributions are rather broad and have 
a considerable overlap. In contrast, the NOMAD ex- 
periment working at significantly higher energies (be- 
tween about 5 and 100 GeV) observed no such excess 
of measured over expected quasielastic cross section [7] . 
Both the MiniBooNE and the K2K experiments could 
obtain good fits to their data in the Fermi-gas model 
only when the axial mass was considerably increased to 
M A = 1-23 GeV [3] or even M A = 1.35 GeV; the latter 
was obtained from a shape-only fit to da/dQ 2 [?]■ The 
sizable increase in the axial mass needed to describe the 
data cannot be ascribed to deficiencies in the Fermi-gas 
model alone. Indeed, Benhar et al. showed that a model 
based on state-of-the-art nucleon spectral functions re- 
quired an even larger axial mass (Ma = 1.6 GeV) for a 
fit of the differential data [5] . 

In contrast to electron scattering, the energy of the 
incoming neutrino is not fixed due to the broad energy 
profile of the neutrino beam; only energy and angle of the 



outgoing lepton can be determined. The cross section at 
fixed energy and scattering angle of the outgoing lepton 
then picks up contributions from different kinematical re- 
gions and thus quite different processes [SJ. For example, 
in the MiniBooNE energy regime there is a strong entan- 
glement of true QE scattering and pion production [10] . 
Events in which a pion produced in the initial interaction 
is absorbed while traveling through the nucleus can also 
lead to knock-out nucleons without any pions and thus 
look very similar to genuine QE events. In Refs. QUI E] 
we have shown that these events can affect the energy 
reconstruction and the extraction of oscillation param- 
eters from neutrino flux comparisons. However, in the 
data mentioned these types of events had already been 
removed with the help of event generators. While this 
removal may not be perfect, the inaccuracies connected 
with it are probably not large enough to explain the ob- 
served excess. 

Here it is interesting to note that also the explicit pion 
cross sections obtained by MiniBooNE are by about a fac- 
tor 1.5 - 2 larger than those obtained with the very same 
event generator (NUANCE) used to analyze the Mini- 
BooNE data @]. A similar result (> 50% excess) holds 
for simulations with the theory-based generator GiBUU. 
In Ref. [T2] 02] we have discussed that the pion yields 
measured by the MiniBooNE experiment Q31 HI] consid- 
erably exceed the pion yields calculated in the impulse 
approximation; the data agree essentially with the calcu- 
lations before any of the strong pion-nucleus final-state 
interactions (FSI) are taken into accountF] The most 
obvious explanation for the common excess in QE scat- 
tering and pion production could be an inaccuracy in the 
flux determination: a 30% higher flux would bring both 
data sets in much better agreement with theoryjj The 
MiniBooNE experiment, however, gives a possible inac- 



1 Qualitatively the same result is found in calculations with the 

neutrino event generator NUANCE [T4llT5| . 

2 The shape of da/dQ 2 could be explained by RPA correlations 
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curacy of only about 10% for the flux [T7], so that this 
explanation seems to be - at least partly - ruled out. 

One major difference of the older data on elementary 
targets is that the mentioned experiments use nuclei as 
target material; in the case of the MiniBooNE experiment 
this is oil (CH2), while for K2K it was water (H2O). It 
is, therefore, tempting to assume that the observed ef- 
fect is due to some nuclear in-medium effect on the ax- 
ial current. Both the analyses of MiniBooNE and K2K 
changed only the axial coupling, leaving the vector cou- 
pling unchanged. In a similar class of models belongs an 
analysis in which the functional form of the axial form 
factor was fitted to the MiniBooNE data, i.e., the dipole 
ansatz for the axial form factor was no longer made. This 
leads to an axial mass M4 = 0.85-1 GeV (defined via the 
derivative of the formfactor at Q 2 = GeV 2 ) [TB]. In 
contrast, a recent explanation of the observed excess by 
Bodek et al. [19] is based on the observation that the 
transverse response of nuclei in inclusive electron scat- 
tering is underestimated by the Fermi-gas model. By 
fitting a nuclear transverse-enhancement factor for the 
electron-data, i.e., by changing only the vector current, 
the authors of Ref . [TH] then are able to describe the Mini- 
BooNE data vs. reconstructed energy and even - approx- 
imately - the disappearance of the excess in the NOMAD 
experiment. It has recently been shown [20] that this ad- 
hoc ansatz of [19: also describes the double-differential 
cross sections. Finally, Meucci et al. [21 find within the 
relativistic Green's function model (RGF) , which is basi- 
cally related to the (one-particle) impulse approximation, 
a very strong effect of FSI. These authors are able to de- 
scribe the MiniBooNE data perfectly well when includ- 
ing the FSI between the ejected nucleon and the residual 
nucleus. In this case the FSI contain implicitly effects of 
other reaction mechanisms. This is similar in spirit to the 
work of Ref. [22] where the authors absorbed higher ex- 
citations beyond lp-lh into the FSI by using Feshbach's 
projection formalism. 

All the explanations described in the preceding para- 



graph rely on the so-called impulse approximation in 
which the incoming neutrino interacts with one nucleon 
at a time only. Early work on electron scattering had also 
shown the importance of random-phase approximation 
(RPA) correlations in the QE-peak region, where these 
correlations tend to lower the cross section [23]; Kim et 
al. showed that similar effects also appear in neutrino- 
induced reactions (24) . Furthermore, from electron scat- 
tering we know that also more complicated processes take 
place, in which the incoming electron interacts with two 
nucleons at the same time. Indeed, in inclusive inelastic 
electron scattering on nuclei these so-called two-particle- 
two-hole (2p-2h) processes become significant at larger 
energy transfers beyond the QE peak, in the so-called 
dip region between the QE peak and the A peak and 
under the A resonance where the strength could not be 
explained in a one-particle picture alone. A rather com- 
plete summary of early attempts for electron scattering 
along these lines can be found in Chapter 5 of Ref. |25j . 
Later work concentrated on a correct relativistic treat- 
ment of meson exchange currents (MECs) the dip region, 
where these indeed lead to a marked contribution to the 
transverse response [26 . 

As a consequence of the experimental method to iden- 
tify QE events in the MiniBooNE experiment (knock- 
out nucleons are not observed) the measured cross sec- 
tion for neutrinos can indeed also contain contributions 
from 2p-2h excitations. Delorme and Ericson noticed al- 
ready in the context of old bubble chamber experiments 
that two-particle-two-hole excitations could contribute 
to the total nuclear response [27]. Following this sug- 
gestion Marteau [28] , in analyzing earlier neutrino ex- 
periments, included 2p-2h excitations and the RPA in 
his analysis. At about this time also Bleve et al. used 
a similar ansatz in their study of the nuclear response 
to neutrinos in the QE region [22]. Martini et al. [29] 
were the first to realize that also the MiniBooNE experi- 
ment could not separate out the QE process, because the 
experiment is insensitive to any outgoing nucleons and 



that 2p-2h processes could thus contribute to the mea- 
sured quasielastic-like cross section. By combining the 
RPA with a calculation of 2p-2h contributions these au- 
thors obtain a good description of the MiniBooNE data 
[291 - I3T] . As expected, the RPA correlations have most ef- 
fect at forward angles where the squared four-momentum 
transfer Q 2 to the nucleus is small. They die out with 
increasing angle and with decreasing muon energy, i.e., 
increasing energy transfer. The calculations do not de- 
scribe the falloff of the QE excess at the higher NOMAD 
energies, but they make predictions for the antineutrino 
cross sections; for these the 2p-2h contributions are sig- 
nificantly smaller than for the neutrino ones. 

The idea of Martini et al. has been taken up by vari- 
ous authors who try to improve on the detailed theoreti- 
cal ingredients. Nieves et al. have obtained in Ref. 32 a 
very good agreement with the energy-separated data, but 
the QE contribution shown there relies on non-relativistic 
approximations which are not reliable at the higher neu- 
trino energies. This can also be seen by comparing their 
QE cross section a(E v ) with the results of other models 
in Ref. [55] : a corrected result gives still a good, but 
somewhat less perfect agreeement for these data with 
a clear discrepancy towards the lower neutrino energies 
[51] , In another calculation by the same authors that 
employs the correct relativistic corrections, but leaves 
out any final state interactions (FSI) 35 , these authors 
have obtained a very good description of the measured 
inclusive double-differential cross sections, when simulta- 
neously readjusting the experimental flux within its in- 
accuracies. Their calculations build on earlier work on 
electron-induced reactions [36, 37 and provide a consis- 
tent theoretical framework; the corresponding electron 
data, though mostly at somewhat lower energies than rel- 
evant for the neutrino experiments, are described quite 
well. As in the work of Martini et al. the disappear- 
ance of the excess when going up to the NOMAD energy 
regime cannot be studied. For antineutrinos these au- 
thors predict, contrary to the results of Martini et al., a 



sizable effect of 2p-2h excitations so that the antineutrino 
cross section becomes nearly as large as that for neutri- 
nos. The FSI effects, which are neglected in Ref. [35] . are 
estimated to be small (< 7%) for the total cross section, 
but could be larger for the lower neutrino energies. The 
fact that Meucci et al. [21] obtain a good description of 
the energy-separated MiniBooNE data due to strong FSI 
within an RGF calculation just underlines the fact that 
different theoretical ingredients may lead to the same fi- 
nal agreement with the flux-averaged neutrino data. 

Finally, an interesting approach is that by Amaro et 
al. [HI [35] which starts with a phenomenological model 
for the neutrino interactions with nuclei that is based 
on the superscaling (SUSA) behavior of electron scat- 
tering data. The scaling function thus determined is 
then directly taken over to neutrino interactions with 
some 2p-2h contributions due to meson exchange cur- 
rents added; those are so far neither gauge invariant nor 
do they contain the axial contributions. There are no 
RPA correlations explicitly taken into account, but they 
may be contained in the scaling function. The overall 
effect of using this model is an increase of the inclusive 
cross section for neutrinos; at forward muon angles the 
calculations come close to the data, but the MEC con- 
tributions die out fast with increasing angle so that the 
cross section is significantly underestimated at backward 
angles. As a consequence the energy-separated cross sec- 
tion obtained for the MiniBooNE experiment - while be- 
ing higher than that obtained from SUSA alone - still 
underestimates the experimental result even when 2p-2h 
contributions are added. Recently, a strong difference be- 
tween neutrino and antineutrino cross sections has been 
obtained within this model, with the 2p-2h effects being 
significantly larger for antineutrinos than for neutrinos 
|40) . The latter result is in contradiction to that ob- 
tained by Martini et al. and by Nieves et al. For short, 
but rather comprehensive discussions of 2p-2h effects see 
alsoRefs. 02112]. 

In summary, there exist various, mutually exclusive ex- 



planations of the observed QE excess. The lp-lh models 
disagree in attributing the observed effects either to the 
axial coupling [31 [TS] , the vector coupling [19] , or to final- 
state interactions [21j . Alternatively - or in addition - 
two-body interaction effects, that are present in electron 
scattering most likely also play a role in neutrino exper- 
iments. Various models agree in their importance, but 
still suffer from problems connected with their treatment 
of gauge invariance, final-state interactions, uncertainties 
in relativistic corrections [311 1351 139] , and upper energy 
limits for their applicability [35 . As a consequence, the 
precise strength of 2p-2h processes and their energy de- 
pendence in neutrino experiments is still uncertain. 

We think that progress can be made by following either 
one of two different paths: 

• On the theoretical side, the calculations of 2p- 
2h contributions can be improved so that gauge- 
invariant calculations of nuclear correlations and 
meson exchange currents become available and 
RPA effects and final state interactions are consis- 
tently calculated with the same forces. At the end, 
for any comparison with experiment, the calculated 
cross sections have to be averaged over quite broad 
energy distributions thus possibly wiping out many 
details of the interaction matrix elements. 

• A more phenomenological approach suggests itself 
in which detailed calculations of 2p-2h contribu- 
tions are avoided by starting with a flux-averaged 
matrix element. Such an approach, if incorporated 
into an event generator, can lead faster to never- 
theless reliable results on experimental observables 
beyond the inclusive cross sections. 

While theoretical consistency arguments may favor one 
model over others the ultimate answer can come only 
from experiments. Inclusive measurements are obviously 
not very sensitive to details of the reaction mechanism 
because they always represent sums over various reac- 
tion mechanisms. For neutrinos this is even more so be- 



cause experiments always involve an average over incom- 
ing energies. A clarification of the reaction mechanism 
must, therefore, come from experiments other than just 
inclusive reactions; exclusive or semi-inclusive reactions 
can give more specific information. If one wants to go 
beyond inclusive cross sections both of the approaches 
named above must be combined with an event genera- 
tor that contains a very good, well tested treatment of 
FSI. In addition, the data necessarily always involve also 
some pion degrees of freedom, even if no final-state pions 
are observed [TU]. Because of the broad energy distri- 
butions in the neutrino beams one thus has to develop 
and use methods to describe consistently scattering pro- 
cesses at low and high momentum transfers with different 
reaction mechanisms (lp-lh, 2p-2h, QE, resonance exci- 
tation, and deep inelastic scattering (DIS)) because all of 
these can contribute to the measured cross sections |43j . 

In this paper we follow the second path from the two 
named above. For the main question raised here (lp- 
lh vs 2p-2h) it is obvious to try to explore possible ob- 
servable consequences of the presence of two-body inter- 
actions in the numbers and spectra of knock-out nucle- 
ons. We will, therefore, investigate the impact of 2p-2h 
processes on knock-out nucleons in a simple model that 
avoids the quite significant difficulties of calculating the 
MEC and correlation effects by making a physically moti- 
vated ansatz for the relevant matrix element. The model 
aims for a first exploration of what will happen if two- 
body interactions play a noticeable role; in a later stage 
we could then implement the results of sophisticated cal- 
culations. We also analyze the impact of 2p-2h processes 
on the reconstruction of the neutrino energy. 

While in this paper we deal only with charged current 
(CC) events it is obvious that such studies are most rel- 
evant for neutral current (NC) interactions where only 
outgoing hadrons are available for the energy and Q 2 re- 
construction. 



II. MODEL FOR 2P-2H CONTRIBUTIONS 

A. Basic properties of GiBUU 

The starting point for our studies is the GiBUU 
(Giessen Boltzmann-Uehling-Uhlenbeck) model that has 
been developed originally as an event generator for 
heavy-ion reactions and then extended to descriptions of 
reactions involving elementary projectiles impinging on 
nuclei. For all the various applications, the very same 
physics input and code are used; this sets the GiBUU 
model apart from many other event generators. A com- 
plete description of the GiBUU model can be found in 
a recent review |44] which also contains a complete list 
of references to applications of the GiBUU model and a 
discussion of results obtained with it. 

In the GiBUU model the spectral one-particle phase- 
space distributions, F(x,p), of all particles are obtained 
by solving the coupled Kadanoff-Baym equations @S] for 
each particle species in their gradient-expanded form [46] 

VF(x,p)-tr{Tf,nS rct (x,p)} vh = C(x,p) , (1) 



with 



VF = { Po -H,F} 



P b 



(2) 



Here {. . .} p b denotes a Poisson bracket. In the so-called 
backflow term [second term on the left-hand side in 
Eq. ffl], which is essential for off-shell transport, f(x,p) 
is the phase-space density related to F by 



F(x,p) = 2wgf(x,p)A(x,p) 



(3) 



where A(x,p) is the spectral function of the particldjand 
g is the spin-degeneracy factor. For on-shell particles 
[A = 5(p° — E)] the phase-space density / is connected 
to the nuclear density by 

d 3 p 



p(x) = g 



(27T)-- 



-f(x,p) ■ 



(4) 



The quantity T in the backflow term is the width of 
the spectral function, and S TOt (x,p) denotes the retarded 
Green's function. Off-shell transport is thus included 
and leads to the correct asymptotic spectral functions 
of particles when they leave the nucleus. The expression 
C(x,p) on the right-hand side of Eq. (fTl) denotes the col- 
lision term that couples all particle species; it contains 
both a gain and a loss term. For a short derivation of 
this transport equation and further details we refer the 
reader to Ref. [4"4"] . 

Because the GiBUU model was developed to be used 
for many different reaction mechanisms, it provides a 
consistent framework for the description of broad-energy 
beam neutrino experiments. The GiBUU model has 
undergone extensive testing with various reaction types 
(see, e.g., Refs. [3H E7]); for the purpose here its val- 
idation with photo-produced mesons on nuclei is most 
relevant. In Refs. [48l|49] we have presented the first ap- 
plications of the model to neutrino reactions with nuclei 
and have analyzed both QE scattering and pion produc- 
tion. The calculations include the FSI of all particles; fur- 
thermore we use relativistic kinematics throughout. The 
procedure for QE and pion-production events is as de- 
scribed in Ref. [48;; further details on the GiBUU model 
can be found in a recent review [44] , Relevant for the 
following discussions is that the QE cross section is cal- 
culated with an axial mass of Ma = 1 GeV. At the 
energies relevant for MiniBooNE, pion production pro- 
ceeds overwhelmingly through the A resonance; its cou- 
pling strength is given by partial conservation of the ax- 
ial current (PCAC) and its form factor is described by a 
modified dipole form factor, as explained in Refs. [48" l l49 | . 



A is normalized as L°° A(x,p)dpo = 1. 



For the MiniBooNE data discussed here, the model, 
which was so far based on the impulse approximation 
and did not contain any two-body interactions for the 
incoming neutrino, also - as all the other models |33j - 
underestimates the measured QE cross sections. We have 
now extended this model by including such processes so 
that now a generator is available that contains the 2p- 



2h effects. The extension will be described in the next 
subsection. 



B. Inclusion of two-nucleon interactions 

In Ref. [H] we have described how we treat the CC 
neutrino interactions with nuclei within the impulse ap- 
proximation. Here we now develop the relevant expres- 
sions for an extension to 2p-2h processes. The starting 
point is the triple-differential cross section for the reac- 
tion v(k) + A -)• l~(k') + X 



dV = |k'| G 2 2 
dfi' dE' E v 4tt 2 ' ' 



(5) 



where k = (E v ,]&) and k' = (U',k') are the incoming 
neutrino and outgoing lepton momenta, respectively, SI' 
is the scattering angle of the outgoing lepton and G is the 
Fermi constant. The squared invariant amplitude \T\ 2 — 
Lu\W^ is given by the contraction of the leptonic tensor 
L and the nuclear hadronic tensor W . 

The total cross section for the interaction of a neutrino 
with a nucleus can be related to the collision rate T of 



the neutrino and is given by [50] 
d-V 



Otot 



m^ = wJ Tix ' k)d3x 



(6) 



Here the local density approximation has been used. The 
collision rate m is directly related to the imaginary part 
of the neutrino self/ energy. 

In transport theory the collision rate can be obtained 
from the loss term in the collision term in Eq. (nj). This 
term contains one-body (resonance decay, for example), 
two-body (QE scattering, for example), and three-body 
(2p-2h) terms. For our purpose here only the two- 
body and the three-body terms are relevant so that T 
is the collision rate for two-body and three-body pro- 
cesses r = r' 2 ) + r^ 3 '. The spatial integration in Eq. M 
extends over the nuclear volume. The two-body interac- 
tions contained in T^ 2 ' have been discussed in detail in 
Refs. [48, 49 and are already implemented in the GiBUU 
model. They contain QE scattering and pion production 
through resonance excitation. 

For the three-body collision processes (y + NiN% — > 
V + N3N4) of interest here, the collision rate, which rep- 
resents the imaginary part of the self energy in local den- 
sity approximation, is given by (for details and notation 
see Sec. 3.3 in Ref. HU) 



i(3) 



IT(z,fc) 



F u {x,k) 
1 SiiSza 



(7) 



d 4 Pi 



d 4 P2 



2E V gm9l , J (27r) 4 2p»y (2^) 4 2p! 
xF 1 (x, Pl )F 2 (x, P2 )(2TT) i S^(k- 



d 4 fc' 



d 4 p 3 



d 4 p 4 



■Pi 



(2vr) 4 2A:' 7 (27r) 4 2p° J (27r) 4 2p° 
P2 -k 1 -p 3 -pi) \M v i2^i>3i\ 2 Fi> (x, k')F 3 (x, p 3 )F 4 (x, p 4 ) 



Here |.M„i2-w'34| 2 is the in-medium matrix element, 
squared and averaged over the spin states of the initial 
particles and summed over those of the final particles. £j 
The Fi = F(xi,pi) are the spectral phase-space densities 
of the two nucleons with which the interaction of the W 



I 

boson takes place (for % = 1, 2) or of the outgoing nu- 
cleons (i = 3,4), and F(x,p) = 2irgA(x, p)[l — f(x,p)] 
contains the Pauli blocking. Since we consider a sta- 
tionary nuclear target all these distributions do not de- 
pend on time. The four-momenta of the nucleons are 
given by pi = (p°,Pi) and their energies by Ei(pi) = 



The spinors in Jv[ are normalized according to u'm= 2E 



y^p? +M 2 + U(x,p l ) (i = 1,2,3,4), where M is the nu- 
cleon mass and U is a space- and momentum-dependent 



nuclear mean-field potential. The gi are spin-degeneracy 
factors and the Sij are symmetry factors (Sij — 1/2 for 
pp or nn pairs, 1 for pn pairs). There is no Pauli blocking 
for the final lepton so that /;/ = 0. Here we neglect all 
energies and momenta in the recoil nucleus. 



such as 



The S^ function in Eq. (M) limits the degrees of free- 
dom. Since we have three outgoing on-shell particles 
with in total nine vector components and four energy- 
momentum conserving constraints, only five degrees of 
freedom are left. One could thus evaluate cross sections 



dV 



(8) 



d 3 k'dn 3 

while all the other kinematical quantities, in particular 
the energies of the two knock-out nucleons, are restricted 
by energy-momentum conservation. In theoretical stud- 
ies of photon- and electron-induced two-nucleon knock- 
out reactions from nuclei, even higher-differential cross 
sections were evaluated [5TT - f55] . For the special case of 
zero energy and momentum of the residual nucleus they 
reduce to those calculated here. 

In order to separate the initial states from the final 
states we now define the collision rate, r^jVi of a neutrino 
with two nucleons with momenta p\ and p 2 at point x 
leading to an outgoing lepton with momentum k', 



T NN (k,pi,p 2 ,k') 



1 S 



:S4 



d 4 p 3 



2Eu 9394 J (27r) 4 2p° 

x (2ir) 4 6^{k+ Pl +p 2 



d 4 P4 

(2n)^pl 

-k' -pa- 



Pi) \M V 



12-W34 



1 F 3 (x,p 3 )F4(x,P4) , 



(9) 



so that we can rewrite Eq. (frl) into 

r(3)c u\ c f d 4 ^ 1 f d4 P 2 



= s 



12 



dV 



(27r)4 2p o 

d 4 P2 



- 5 F 1 (x,p 1 )F 2 (x,p 2 ) 



d 4 fc' 1 — 

T NN (k,p 1 ,p 2 ,k') — F v (x,k') 
9i' 



Fi(x,pi)F 2 (x,p 2 ) 



(2-n) 4 2k' 
d 3 fc' 



-T NN (k,p 1 ,p 2 ,k') 



(10) 



(27r) 4 2p° J (2tt) 4 2i>S* "">*"■'* ^->*"» J (2n) 3 2E' 
In going from the first to the the second line we have used for the outgoing lepton Fy = 2wgi/d(k' — E') so that in 



Eq. (10) the integration over k' can be carried out, reducing the integral over k! to a three-dimensional one. The 



differential cross section can now be obtained from Eqs. ffity and (10 1, 



d 3 CT (3) 



E„ 



1 



rSi 



A A x 



d 4 pi 



d 4 p 2 



o F 1 (x,p 1 )F 2 (x,p 2 )T NN (k,p 1 ,p 2 ,k') . 



(11) 



d 3 /fc' |k|£'2(27r) 3 ~ lz J " "J (27r) 4 2p? J (27r) 4 2p° 

In this form the actual two-body interaction contained in Tnn is separated from the phase-space distributions of the 
initial particles. This corresponds to the so-called factorization already used by Gottfried [5Tj . 

Up to this point the formalism is quite general as far as the outgoing particles are concerned. These could be either 
NN, NA, or AA pairs. We now simplify the two-body collision rate Tnn- Since we are primarily interested in 
nucleons being ejected into the continuum the Pauli-blocking factors 1 — f 3 A can be neglectecrl so that F = 2ngA. 



In the actual calculations we keep the Pauli-blocking factors, 



thus allowing also for a possible recapturing of the knock-out 
nucleons. 



We now eliminate the four-dimensional momentum-conserving d function in T^jv by integrating over d 4 p 4 . With 
fixed P = k — k' + pi + P2 — p 3 + P4 the function r^jy (f9J) is then given by 



^NN(k,p 1 ,p 2 ,k') = ^— ^<S 34 



d 4 p 3 



f 



_ , ^^^^v^A^x.p^A^x, P - p 3 ) 



(12) 



(27r)42p§2( J P"-pg) 1 

The remaining integral over d p 3 can be further simplified by using the so-called quasiparticle approximation, i.e., 
neglecting the width in the spectral functions of the outgoing particles, i.e., restricting ourselves to outgoing nucleons 
by setting 



Mx,p) = 6(p? - Ei) 



with Ei = yp~f + M 2 . After integrating Eq. ( 12 ) over dp® one obtains 

d 3 p 3 



T NN (k,p 1 ,p 2 ,k') 



c 

E v 3 J (2tt) 3 2E 3 2E 1 



^6(P - E 3 -E±)\M vl2 - 



■Z'34| 



(13) 



(14) 



with E 3 = \ZpI + M 2 and E 4 = \/(P — P3) 2 + M 2 = Pq — E 3 . Integrating out now |p 3 |, exploiting the 6 function, 
gives 

1 



T NN (k,p 1 ,p 2 ,k') 



J_ n r d» 3 p 3 
E u 3 UJ (2n)3E 3 E 4 



d(E 3 +E 4 ) 
d|p 3 | 



\M U 



12->Z'34| 



eZ Sm ^ 



pI 



d^3 YJ-. p 

4t ' ,/12 ^' 341 (P - £fe)|p 3 | +S3IP-P3I 



(15) 



Here P3 = | JD3 1 is determined as the solution of the equation Pq — E 3 — E4 = and depends on Po, P and cos(Z(P, P3)). 

axim 

d 3 p 2 



Inserting Eq. (15) now into Eq. (11) gives, after using the quasiparticle approximation also for the initial nucleons, 

d 3 Pl 



dV 3 ) 



dO' dE' 



|k'| c «, 1 

M 919 3 * "(2^8 
dft 



d 3 x 



(2 7 r) 3 2 J Bi 7 (2n) 3 2E 2 



fi{x,pi)f2(x,p 2 ) 



pi 



^ lM ^ l ' 3il2 IP^E 3 )\p 3 \ + E 3 \P-p 3 \ 



(16) 



where p° = Ei in /j. After inserting the spin-degeneracy factors gi = g 2 = 2, using S34S12 = 1/2 for CC reactions 
(the initial pp state is not possible for CC reactions) and assuming isospin independence of the CC interaction matrix 
element and the phase-space distributions, we obtain finally for the triple-differential cross section 

dV 3 ) ik' 1 1 r ,» f d 3 Pl r d 3 P2 



dWdE' |k| 4(2tt) 4 

d^ 3 



d 3 x 



(27r) 3 2^i J (2ir) 3 2E 2 



fl(x,Pl)f2(x,p 2 ) 



it 



An |M, 12 ^34| 2 (po _^ 3)| . 3 +rj |p_ p! 



(17) 



Higher-fold differential cross sections such as Eq. p| can be obtained similarly. 



Comparison of Eq. (17) with Eq. (|5|) allows one to express the nuclear tensor W^ x for the nuclear systems in terms 



of that for the 2N system, 



W» = Sufd*z\[^f^Mz tPl )M*,P2)Fw* 



(18) 



(2tt) 3 J (2tt) 3 

Here F = (pnn ■ k)/{p NN E v ) is the usual flux factor that transforms the cross section from the two-nucleon (NN) 
to the nuclear system, and w^ x represents the hadronic tensor for the weak interaction with the two-nucleon system. 
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We neglect the Lorentz transformation from the moving NN system to the stationary nuclear system. The matrix 



element |.Mi/i2->-z'34| 2 and the hadronic tensor w^ are related by 

|.M„12->1'34| 2 = L^w» x , 



(19) 



where L^\ is the leptonic tensor. The phase-space distributions / appearing here are solutions of the transport 
equation (fTl) and thus contain the effects of the nuclear potentials as well as of all final-state interactions. Starting 
from a reliable model for neutrino interaction with the 2N system one could calculate w^ x consistently in a Fermi-gas 
model. 



Eq. (17| shows nicely [cf. Eq. Q] that the 2p-2h pro- 
duction cross section is oc p 2 where p is the nuclear den- 
sity. The fact that both initial nucleons are at the same 
location x reflects spatial correlations between them. The 
dependence on the p 2 also immediately implies that FSI 
will be very important for this process, both for the inclu- 
sive cross section and the semi-inclusive one for knock-out 
nucleons. 



III. RESULTS 



There are various detailed calculations for the ma- 



trix element M. appearing in Eq. (17) presently going 
on jUmmmil]. These calculations are quite de- 
manding and still have to fight with problems such as 
gauge invariance, relativity, the inclusion of final state 
interactions, and limited energy regions of applicability. 
For example, the authors of Ref. [35] give as an upper 
limit of validity of their calculations an energy of 1.0 - 
1.2 GeV, whereas the MiniBooNE flux extends well be- 
yond this limit. Since the main purpose of this paper is to 
look for possible distinctive experimental observables for 
the presence of 2p-2h processes and to investigate their 
influence on the energy reconstruction, we now make sim- 
plifying assumptions that, as the results will show, still 
seem to capture most of the relevant physics. The start- 
ing point for our approximation is the observation that 
naturally all neutrino cross sections with nuclei are av- 
eraged over the incoming neutrino energy distribution. 
For the double-differential cross section d 3 a/dil'dE' at a 



fixed E' this directly translates into an average over en- 
ergy transfer and Q 2 that will wipe out many details of 
the matrix element. We, therefore, directly parametrize 
an average matrix element in terms of the two Mandel- 
stam variables s and t = —Q 2 and fit the s dependence 
to the energy-separated MiniBooNE data in order to fix 
the overall strength of this contribution at each energy. 
With this input we then calculate, as a first test, the 
double-differential cross sections and compare them with 
the MiniBooNE data. Furthermore, we investigate the 
numbers and spectra of knock-out nucleons. Finally, we 
will discuss the implications for the energy reconstruc- 
tion. 

We have worked with various parametrizations and will 
give results here only for two. One is a pure phase-space 
model, in which A4 is constant (model I); the constant 
is fitted to the difference of the energy-separated Mini- 
BooNE data and the calculated QE cross section. In the 
other model we approximate the hadronic tensor w by 
the transverse projector Pj — — g M „ — q^q v /Q 2 (model 
II), again fitting the overall strength and an s depen- 
dence, and then contract it with the lepton tensor. The 
latter model is meant to mimic the prevalence of two- 
body reaction components which are known to be trans- 
verse [5B]. In all these studies RPA effects are not taken 
into account. Since these tend to lower the cross sec- 
tion in particular at forward angles, our calculations may 
somewhat underestimate the contributions of 2p-2h ef- 
fects there. We note, however, that very recent studies 
show that RPA effects are neglible for neutrino energies 
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larger than 1 GeV and are much smaller than theoretical 
uncertainties already for E v > 0.7 GeV [57] . 

All calculations are done for charged current interac- 
tions on a 12 C target, using an axial mass Ma = 1 GeV 



for QE processes. We evaluate Eq. (17 1 by picking ran- 



domly a pair of nucleons (1,2) at the same location with 
momenta chosen out of the Fermi sea; all further correla- 
tions are assumed to be contained in the matrix element 
A4 . Then — at fixed outgoing lepton energy and angle — 
we choose two outgoing nucleons (3,4) that obey energy 
and momentum conservation, i.e., that are weighted ac- 



cording to the integrand of the H3 integration in Eq. ( 17 ) 



We then follow these two nucleons through and out of the 
nucleus with all the final-state interactions (elastic and 
inelastic scattering, charge transfer) included. In doing 
so we assume that the density and potential of the tar- 
get nucleus are not significantly disturbed (the so-called 
frozen approximation). The whole event finally receives 



a weight given by the cross sections (17). This procedure 



yields both the knock-out cross sections and — after in- 
tegration over the momenta of the outgoing nucleons — 
also the inclusive double-differential cross section. We 
note that the calculations include the FSI of all particles 
and use fully relativistic kinematics throughout. 

A. Inclusive cross sections 

In Fig. [l] we show for reference again the MiniBooNE 
data 0] (those with the removed contribution from the 
absorbed pions) together with the predictions of the 
GiBUU model for the QE cross section [48,. The dif- 
ference between the data and our QE calculation gives 
the 2p-2h contribution which is fitted here by model I. 
It is obvious that this very simple phase-space model is 
not perfect. It gives a cross section that rises too steeply 
at energies above about 1 GeV (but is still within the 
experimental errors bars). 

Fig. [2] shows the calculated cross section da/dQ 2 in 
model I. While the gross structure is quite well repro- 
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FIG. 1. (Color online) MiniBooNE data for QE scattering 
(from Ref . [4] ) together with the result of the GiBUU model 
for the integrated QE cross section (dash-dotted line), calcu- 
lated with Ma = 1 GeV. The dashed line gives the contri- 
bution of 2p-2h excitations using a constant matrix element 
(model I, see text); the solid line gives the sum of both. All 
cross sections are per neutron. 

duced, there is a clear disagreement at lower Q 2 . Taking 
into account the RPA correlations would cure this prob- 
lem. Indeed, as it was shown in Ref. |31) . RPA correla- 
tions suppress the true-QE cross section by 20%-25% at 
Q 2 < 0.2 GeV and by a smaller percentage further up to 
Q 2 < 0.4 GeV. We show this result here only to exhibit 
the Q 2 dependence of the 2p-2h contribution on phase 
space alone. Phase space alone thus favors an increase of 
the 2p-2h cross section towards smaller Q 2 ; the sudden 
fall-off towards Q 2 — is caused by Pauli-blocking and 
phase-space restrictions due to the finite muon mass. 

In order to understand this result we show in Fig. [3] 
the various contributions to the double-differential cross 
section for four different angles. The cross sections are 
averaged over the MiniBooNE flux. It is seen, first, 
that the 'bare' QE contribution describes the forward 
data very well, but comes out somewhat too low at the 
higher angles. It is furthermore seen that the shape 
and over-all size of the 2p-2h contribution is rather in- 
dependent of angle (within a factor of 2), amounting to 
w 0.2 • 10~ 38 cm 2 /GeV in its maximum at T M « 0.2 GeV. 
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FIG. 2. (Color online) MiniBooNE data for QE scattering 
(from Ref . [4] ) together with the result of the GiBUU model 
for the da/dQ 2 cross section per neutron for model I. All 
contributions are averaged over the MiniBooNE flux. 



The overall agreement reached by fitting one number, the 



size of the squared matrix element in Eq. (17), is already 



quite good and certainly better than in much more so- 
phisticated model calculations of the matrix element |38j . 
The main discrepancies appear at very forward angles 
cos 6* = 0.95, where the cross section is overestimated 
due to the absence of RPA in our calculations, and at 
slightly larger angles where the peak cross section is un- 
derestimated. 

Since these results are averaged over the MiniBooNE 
flux it is instructive to look at the inclusive cross sec- 
tion also for a fixed energy calculation. This is shown 
in Fig. HI which gives the cross section for the various 
components at the fixed neutrino energy of 1 GeV as a 
function of energy transfer v. For convenience the in- 
variant mass W, as it would be determined experimen- 
tally from lepton kinematics W 2 = M 2 + 2Mv — Q 2 with 
Q 2 = 2E y {E' —\k'\ cos 9'), is also labeled at the top hor- 
izontal axis of each panel. Fig. HJ shows that the QE con- 
tribution drops drastically when going from a forward to 
a backward angle while the 2p-2h contribution decreases 
only slightly. This results in a significant rise of the ratio 
of 2p-2h/QE with angle; at backward angles both contri- 
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FIG. 3. (color online) The double-differential cross section 
da/(dcosO t j,dT ti ) per neutron for model I for four different 
muon scattering angles versus muon kinetic energy. 

butions become comparable. 

In order to exhibit the influence of a particular Q 2 
dependence of the matrix element we now model the 
hadronic currrent by the transverse projector (model II), 
again fixing the overall magnitude of the matrix element 
by fitting the MiniBooNE energy-separated data. In this 
fit we also allow for an s-dependence of the matrix el- 
ement, which leads to a much better correspondence to 
the data, as shown in Fig. [5] 

In Fig. Kil we show again the results for da/dQ 2 . The 
agreement is now significantly better, with still a little 
too much strength at small Q 2 . However, we note that 
our calculations do not contain any RPA correlations that 
tend to lower the cross section at forward angles; the 2p- 
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FIG. 5. (Color online) MiniBooNE data for quasielastic scat- 
tering (from A.) together with the result of GiBUU for the 
total cross section per neutron. As Fig. [I] but now for model 
II: the hadronic tensor has been assumed to be proportional 
to the transverse projector and contains an s dependence, see 
text. 
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FIG. 4. (Color online) The cross section dcr/(d^dcosS / j) per 
neutron for model I for four different angles as a function of 
energy transfer v. The incoming neutrino energy is fixed at 
lGeV. 



FIG. 6. (Color online) MiniBooNE data for quasielastic scat- 
tering (from Ref. [4j) together with the result of the GiBUU 
model for the da/dQ 2 cross section per neutron for model II 
(see text). 



2h strength should thus probably be even a little larger 
at forward angles than in our studies here. 

In Fig.[7]we show the double-differential cross sections 
for all angular bins, as obtained in model II. Immediately 
evident is the fact that the QE process alone describes 
the forward-angle cross sections quite well, but becomes 
significantly too low at backwards angles. Conversely, 



the 2p-2h contribution is seen to be very small at for- 
ward angles and grows with angle so that at backwards 
angles it is as large as that of the QE process. It is 
seen that at forward angles the 2p-2h contributions are 
largest at small T M , i.e., at large energy losses below the 
QE peak, whereas at larger angles the 2p-2h energy loss 
peaks roughly at the position of the QE peak. Overall, 
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the description reached in our simple model using the 
MiniBooNE flux is obviously quite good and certainly 
comparable to that obtained in a microscopic model with 
flux renormalization (cf. Fig. 1 in Ref. [55]). 

This agreement with the results of Ref. [35] as well 
as the comparison of the results for model I (extreme, 
phase space only) with those for model II (some physics 
basis, transverse) shows that the flux-averaging leads 
to a remarkable insensitivity of the inclusive double- 
differential cross sections to the input matrix element. 
In general good agreement is reached with these simple 
parametrizations of the interaction matrix elements. It 
thus seems to be obvious that by a proper parametriza- 
tion of the (s, Q 2 ) dependence of the matrix element a 
nearly perfect agreement for all the observable inclusive 
cross section data could be reached. While these inclu- 
sive double-differential cross sections are 'cleanest' in the 
sense that they contain only a small model dependence 
(in the removal of pion production events in which the 
pion is absorbed in the target) they also seem to be fairly 
insensitive to details of the hadronic tensor. In particular 
they do not allow one to separate QE, RPA, and 2p-2h 
effects in a unique way. 

We abstain here from trying to get a perfect fit for 
the matrix elements and instead now look for observable 
consequences of the 2p-2h processes, that go beyond just 
inclusive cross sections with their inherent ambiguities as 
far as their detailed composition is concerned. 

B. Nucleon knock-out 

In this section we thus discuss now properties of knock- 
out nucleons and their sensitivity to the underlying mech- 
anism (lp-lh vs 2p-2h). Hereinafter the calculations 
contain all relevant processes, i.e., QE, resonance and 
background pion production, and model II for 2p-2h pro- 
cesses. 

In Fig. [8] we show the knock-out cross sections for 1, 
2, 3, and 4 nucleons before (dot-dashed curves) and after 



(thin solid curves) FSI. Various original processes (2p- 
2h, QE, A, higher resonances, 1-pion background) can 
contribute to the same multiplicity. For the cross sec- 
tion after FSI some of them (2p-2h, QE, A) are explic- 
itly shown. There is no interference between different 
contributions to the same multiplicity event because the 
reactions are all semi-inclusive. 

Before FSI, 1-nucleon knock-out comes from all the 
processes considered in the model. This includes a tiny 
contribution from 2p-2h processes, which may happen 
when one of the outgoing nucleons has a too-small kinetic 
energy and remains bound inside the nucleus. For 2- 
nucleon knock-out before FSI the contribution is solely 
from the 2p-2h processes; 3- and 4-nucleon knock-outs are 
not possible (recall that we do not consider DIS processes 
here). 

After FSI, all the curves are noticeably different from 
the corresponding ones before FSI. For 1-nucleon knock- 
out FSI lower the cross section, whereas for all the higher 
multiplicities they increase it. In each case the total cross 
section rises with energy, driven by the QE contribution 
at the smaller energies E < 0.5 GeV and by pion produc- 
tion above that energy. Only the A contribution to the 
latter process is explicitly shown in the figure. After an 
initial rise the 2p-2h contribution stays rather constant. 

For the 2-nucleon knock-out, the QE contribution is 
even slightly bigger than the 2p-2h one. This is a FSI 
effect in which after the initial QE process the outgoing 
nucleon collides with bound nucleons in the nucleus and 
knocks out one of them. At the same time the A con- 
tribution rises and becomes comparable with both at a 
neutrino energy of about 1 GeV. Two processes are re- 
sponsible for this rise. First, the main nonpionic decay 
channel of A in medium, namely radiationless transition 
AJV — » NN, leads to the 2 outgoing nucleons. Second, 
both pion and nucleon produced in A decay can collide 
with bound nucleons in the nucleus and knock out one of 
them. This would again give 2 outgoing nucleons. Those 
again can collide with bound nucleons knocking out the 
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FIG. 7. (Color online) The double-differential cross section do"/(dcos# M dT M ) for various muon scattering angles versus muon 
kinetic energy for model II. The data are taken from Ref. [J. 



third one, thus leading to an avalanche effect. 

The very same processes increase the number of knock- 
out nucleons even more effectively when started by two 
nucleons originating from a 2p-2h vertex. They are re- 
sponsible for the significant decrease of 2p-2h cross sec- 



tion in the 2-nucleon channel and the increase of that in 
channels with higher multiplicity. After FSI, for all 2-, 3- 
and 4-knockout nucleons the relative 2p-2h contribution 
amounts to about 20%. The relative 2p-2h contribution, 
as compared to the QE one, increases with multiplicity, 
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FIG. 8. (Color online) Cross-section for multi-nucleon knock- 
out as a function of neutrino energy. Here the hadronic tensor 
has been assumed to be proportional to the transverse projec- 
tor (model II). In each case the total cross section is given by 
the thin solid line. The QE contribution is indicated by the 
long-dashed line. The short-dashed line shows the A contri- 
bution, while the 2p-2h contributions is indicated by the thick 
solid line. All these lines show the results after FSI, while the 
dot-dashed curves shows the total cross section before FSI. 

which is also a consequence of the higher starting multi- 
plicity (2 nucleons as compared to 1). 

Thus, FSI are strong enough to lead to avalanching 
that masks the starting event. Indeed all the 3- and 4- 



nucleon knock-out events are due to FSI. 

On the other hand, the topology of the final state of 
the first, primary interaction (1 versus 2 outgoing nu- 
cleons) is sufficiently different for lp-lh and 2p-2h pro- 
cesses. Even the avalanchig from FSI does not wash out 
these differences: 2p-2h processes hardly contribute to 
the 1-nucleon knock-out. This implies that the method 
to absorb many-body effects into a rescaling of the mag- 
netic (single-particle) electromagnetic form factor of the 
nucleon leads to incorrect final states if used in a stan- 
dard event generator as suggested by Bodek et al. [T5] . 

There are three interesting messages contained in 
Fig. [8| First, QE scattering and A contribute signifi- 
cantly not only to 1-, but also to 2-, 3- and 4- nucleon 
knock-out; 2p-2h processes noticeably contribute to 2-, 
3- and 4-nucleon knock-out. Thus, particle multiplicities 
alone do not distinguish between initial lp-lh and 2p-2h 
excitations, which means that comparison with theory 
has to be quantitative and not just qualitative. 

Second, after FSI 2p-2h processes contribute less than 
30% to the 2-nucleon knock-out. Thus, even a perfect 
theoretical model for 2p-2h interactions, taken on its own, 
is not sufficient to predict the nucleon output. A realistic 
description of FSI is essential for any verification of the 
reaction process through comparison with experimental 
data. 

The third message is an expected one: since the 1- 
nucleon channel is practically exclusively due to the ini- 
tial QE process and pion production, gating on events 
with one and only one proton (and no pions) can thus 
help to select true QE scattering. We have shown in 
Refs. [TQl [TT] that even then detector thresholds affect 
the measured QE cross section significantly. For exam- 
ple, in a tracking detector with a typical detection thresh- 
old of proton kinetic energy of 0.2 GeV only about one 
half of the full QE cross section is actually measured; the 
rest has to be reconstructed with the help of an event 
generator. 

Finally, in Fig. [9] we show the kinetic energy spectra 
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of knock-out nucleons for events with 1 proton and any 
number of neutrons, 2 protons and any number of neu- 
trons, at least 1 proton and any number of neutrons, and 
2 neutrons and any number of protons in the final state 
for the MiniBooNE flux. First of all, the plots show that 
the importance of 2p-2h knock-outs relative to others in- 
creases dramatically with the number of ejected protons. 
While it is insignificant for the "1 proton X neutrons" 
events it is of comparable size for the higher-multiplicity 
events. Second, there is a pileup of cross section at small 
kinetic energies Tjy (below about 0.1 GeV) due to FSI; 
the 2p-2h contributions alone exhibit a similar behavior, 
though with a slightly smaller rise at small Tjy. The 
2p-2h contributions increase the cross sections at small 
kinetic energies significantly, by a factor 2-3. 

As already discussed, of primary interest for the mod- 
ern experiments are the QE-like events, which by ex- 
perimental definition for Cerenkov detectors comprise 



all events with no pions in the final state. In Fig. 10 



we show the nucleon knock-out cross sections for these 
events. Similar to the previous discussion, wc compare 
the calculations before and after FSI. 

For the final state with 1 proton and neutrons the 
by-far dominant contribution comes from true QE scat- 
tering. A tiny contribution from the 2p-2h processes (as 
discussed above, this is possible if one of the outgoing 
nucleons has a very low kinetic energy) is negligible. FSI 
decrease the cross section; this is because a higher-energy 
nucleon may rescatter and knock out more nucleons from 
the nucleus. For the final state with protons and 1 neu- 
tron the only contribution before FSI comes from the 2p- 
2h processes, if one of the nucleons remains bound. FSI 
of outgoing protons coming from the intial QE process in- 
crease this cross section by a factor of 5. The final state 
with 1 proton and 1 nucleon before FSI is again solely of 
2p-2h origin. FSI noticeably change the weight of var- 
ious contributions: the output from 2p-2h is reduced, 
while QE process contributes significantly. A final state 
with 1 proton and 2 neutrino is only possible after FSI. 
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FIG. 9. (color online) Kinetic energy spectra for knocked-out 
protons (three upper panels) and neutrons (bottom panel) for 
events with any number of mesons in the final state for the 
MiniBooNE flux. 



Here one also sees a mixture of events originating from 
QE, A, and 2p-2h. 

Thus, the only clear signature of true QE scattering 
are the events with pions, 1 proton, and neutrons in 
the final state. The purity of such an event sample would 
be very high, as the nearly negligible difference between 



the dashed and thin solid curves in Fig. 10 shows; how- 
ever, the selection efficiency is only at the level of « 70% 
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C. Energy reconstruction 



FIG. 10. (Color online) Cross section for nucleon knock-out 
as a function of neutrino energy for events with pions in the 
final state. 
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Many neutrino experiments, including MiniBooNE, 
use QE scattering as a tool to determine the incoming 
neutrino energy. This is usually done by using quasifree 
kinematics for the QE process on a neutron at rest. It is 
obvious that Fermi motion is thus neglected from the out- 
set and can be expected to lead to a broadening of the 
reconstructed energy around the true neutrino energy. 
Any binding effects in these analyses are taken into ac- 
count by assuming a constant energy shift. In Ref. [TU] we 
have shown that already the entanglement of QE scatter- 
ing and pion production leads to a significant downward 
shift of the reconstructed energy. In this section we now 
analyze the influence of the initial 2p-2h excitations on 
the energy reconstruction. Based on the implementation 
of 2p-2h effects as described above we have analyzed the 
energy reconstruction on an event-by-event basis. The 
calculations contain all relevant reaction mechanisms, in 
this case QE scattering and pion production, and include 
final-state interactions. 

The relevant formula used for the energy reconstruc- 
tion is based on the assumption of QE scattering on a 
nucleon at rest [17] . 



E v = 



2(M n - E B )E' - (E% - 2M n E B + to 2 + AM 2 



M n - Ef 



Ifc'l cos# L 



(20) 



Here Eb > is a constant binding energy, M n the mass 
of the neutron, and AM 2 = M 2 - M 2 . E', \k'\, and 
9^ are the energy, momentum, and angle of the outgoing 
muon. 

There are two features of this expression that affect the 



analysis of a 2p-2h event using Eq. ( 20 1 , which is correct 



only for true QE scattering, i.e., a lp-lh event. First, in 
the last section we have shown that, at forward angles, 
where the cross section is largest, the 2p-2h contribu- 
tions are largest at small T M , below the QE peak. When 
analyzing such events with the help of the one-body ex- 



pression ( 20 1 this leads to a lower reconstructed neutrino 



energy than the true one. Second, while QE scattering 
is strongly forward peaked, the 2p-2h events are fairly 



flat (within a factor of 2) in lepton angle (see Figs. [3] and 
Ff\. The relatively strong yield at backward angles will 
lead to a larger neutrino energy, in particular for inter- 
mediate muon energies. Since both effects are present we 
expect a fairly flat behavior of the 2p-2h contribution to 
the energy reconstruction. 

In Fig.[TT]wc plot the distribution of the reconstructed 
neutrino energy obtained using the MiniBooNE recon- 
struction method with Eb — 8 MeV (which is a typical 
binding energy in the GiBUU code, as opposed to the 
value 34 MeV used by MiniBooNE) for the fixed true neu- 
trino energies El ruc = 0.5, 0.7, 1.0, and 1.5 GeV. The 
main contribution to all QE-like events is given by the 
true QE events, which shows a prominent peak around 
the true energy. The peak is approximately symmetric 
and has a width of about 0.1 GeV: this broadening is 
caused by Fermi motion. For A— induced events the dis- 
tribution is not symmetric, with a broad peak at lower 
energies. For a more detailed discussion of these effects 
seeRefs. [TO1HT]. 

For all energies the broad 2p-2h contribution is peaked 
around 1/2 of the true energy. For the lower neutrino 
energies of 0.5 and 0.7 GeV most of the total distortion 
is caused by 2p-2h processes because at this low energy 
there is only little A excitation. For the true neutrino 
energies of 1.0 and 1.5 GeV, on the other hand, the con- 
tributions of 2p-2h events are comparably small because 
here A excitation plays a major role and because their 
strength is distributed over the whole energy range from 
0.2 to 1.8 GeV, with a flat maximum again around 1/2 
of the true energy. Both the 2p-2h effects as well as the 
A excitation lead to a shift of the reconstructed energy 
towards smaller values, or, vice versa, for a given recon- 
structed energy the true energy always lies higher than 
the reconstructed one. The effect is most pronounced at 
low true energies and becomes smaller at higher ener- 
gies. Both of these results agree with the recent analyses 
by Martini et al. [SS] and with the very recent results 
of Nieves et al. [57) . Both of these publications do con- 
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tain some treatment of RPA correlations which are most 
prevalent at lower energies E v < 0.7 GeV and at forward 
angles. On the other hand, the work [57] does not con- 
tain any pion (or A) degrees of freedom and thus misses 
part of the lower hump in Fig. |11| This is of no concern 
as soon as calculations are compared to the cross sec- 
tion with subtracted pion-induced QE-like background 
(as done by MiniBooNE), but is essential for any extrac- 
tion of oscillation parameters. 



reconstructed E v 
true E v 
true E v : QE 
true E v : Delta 
true E v : 2p2h 





FIG. 12. (Color online) Event distribution for QE-like pro- 
cesses as denned by the MiniBooNE experiment versus the 
true (upper picture) and reconstructed (lower picture) neu- 
trino energy. The thin solid and dot-dashed curves repeat 
the distributions for true and reconstructed energy in both 
figures. 



This is also seen in Fig. 12 which shows the event dis- 
tribution (also called flux-folded cross section) for QE- 
like processes that would be observed by the MiniBooNE 



experiment. Only part of these zero-pion events are gen- 
uine QE events. The rest can originate from initial A 
production, higher resonance production (not shown in 
the figure), pion-background events (not shown in the 
figure), or 2p-2h processes. Since in the numerical sim- 
ulation we know the true neutrino energy and calculate 
the reconstructed energy as done by MiniBooNE, we can 
predict the event distribution versus both of them. The 
event distribution versus the true energy is shown by the 
curve labeled 'true E v \ Comparing it to the one versus 
the reconstructed energy one sees that there is a system- 
atic distortion in energy reconstruction: the peak of the 
event distribution is shifted now by about 100 MeV to 
lower energies with this shift becoming smaller for en- 
ergies above 1.3 GeV. Qualitatively this same behavior 
was already observed and discussed for the K2K and the 
MiniBooNE flux in Ref. [TU] as a consequence of pionic 
excitations alone. Here, now also the 2p-2h contribu- 
tion is included in this analysis; it leads to the noticeable 
downward shift of the reconstructed energy curve on the 
low-energy side of the maximum where the pion produc- 
tion is still smallrl This distortion is important since it 
directly distorts the oscillation pattern and thus affects 
the extraction of the CP-violating phase Acp • 

So far we have concentrated on a discussion of 2p-2h 
effects in the MiniBooNE experiment, i.e., a Cerenkov 
counter experiment. For completeness we now discuss 
also the K2K experiment which uses a tracking detector, 
which is sensitive to protons but not neutronsjj In this 
experiment an event is identified as QE-like if no pions 
and exactly 1 proton is found in the final state. Evi- 
dently, the proton can be accompanied by any number 



6 We note that there is a slight inconsistency in our analysis be- 
cause we have fitted the matrix element to the cross section vs. 
the reconstructed energy in Figs, [l] and [5] The error connected 
with this is still well within the general uncertainties in the fig- 
ures mentioned. 

7 The presently running T2K experiment employs both techniques 

in its (different) near- and far-side detectors. 
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of unobserved neutrons. As we have shown in Ref. |10j . 
in such a detector there is nearly no distortion from A 
excitation, but the total QE rate is underestimated (be- 
cause of charge exchange in the FSI) and has to be re- 
constructed by means of event generators. 

In this experiment the neutrino energy is reconstructed 



according to Eq. ( 20 1 , but neglecting the binding energy, 



i.e., with Eb = 0. Our results for the cross sections ver- 



sus the reconstructed energy are shown in Fig. 13 One 
can easily conclude, that the experimental possibility to 
detect a proton in the final state and to select events 
with one proton only significantly reduce the contribu- 
tion from the initial A resonance production, while the 
distortion from the 2p-2h processes (downward shift) re- 
mains on the same level as for the MiniBooNE Cerenkov 
detector. 



IV. SUMMARY 

In summary, understanding the reaction mechanism of 
neutrinos with nuclei is mandatory for a reliable determi- 
nation of the neutrino energy. The large excess of QE-like 
events observed by MiniBooNE over those calculated by 
commonly used event generators with the world-average 
value for the axial mass of around 1.0 GeV has there- 
fore led to a multitude of partly contradictory models 
that all aim to describe the same data, starting from 
quite different initial lp-lh or 2p-2h interaction mech- 
anisms. Since broad energy-band neutrino experiments 
necessarily contain an averaging over quite different re- 
action mechanisms, the inclusive data alone do not allow 
for an unambiguous experimental verification of any par- 
ticular model. 

We have, therefore, investigated the effects of 2p-2h 
interactions by extending the GiBUU event generator, 
which contains all the relevant reaction mechanisms, on 
both inclusive and more exclusive observables. Because 
of the inherent kinematical averaging in any broad-band 
neutrino experiment we have approximated the matrix el- 
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FIG. 13. (Color online) Distribution of the reconstructed neu- 
trino energy for QE-like processes as defined by the K2K ex- 
periment (0 pions, 1 proton, and any number of neutrons in 
the final states) versus the reconstructed neutrino energy for 
the true energy of the incoming neutrino of 0.5, 0.7, 1.0, and 
1.5 GeV. 
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ement for neutrino interactions with two nucleons by an 
ad-hoc average ansatz, either using just a simple phase- 
space model or a model in which the hadronic tensor is 
modeled by the transverse projector for the W bosons. 
In both cases quite reasonable descriptions of the inclu- 
sive double-differential cross section depending on muon 
observables only are obtained. We interpret this result 
as an indication that inclusive neutrino cross sections in 
broad energy band experiments are not sensitive to de- 
tails of the vNN interaction as long as average properties 
of the interaction are described well enough. 

Within this same model we have then calculated the 
numbers and spectra of knock-out nucleons. We find that 
indeed 2p-2h primary interactions increase the cross sec- 
tions for multi-nucleon knock-out. Both the numbers 
and the spectra of knock-out nucleons are changed so 
that detailed measurements should be able to pin down 
the strength of the 2p-2h processes. However, even in 
the presence of more complicated initial interactions the 
one-proton knockout channel is only slightly affected by 
2p-2h contributions. Experiments that can actually re- 
strict the number of knock-out nucleons to one and only 
one proton should, therefore, be able to provide a clean 
(with purity at the level of 90%) sample of QE events. We 
have shown earlier that such events also are not signifi- 
cantly contaminated by A excitations 10] ; they do miss, 
however, part of the QE cross section because of charge 
changing FSI on the outgoing nucleons. The selection 
efficiency is expected to be at the level of 70%. 

Finally, we have shown that the energy reconstruction 
is affected by the presence of 2p-2h interactions mainly 
at the lower neutrino energies, while the A excitations 
become more essential at the higher energies. 2p-2h in- 
teractions lead - as the A excitations do - to a down- 
ward shift of the reconstructed energy. Their contribu- 
tion to the reconstructed energy distribution is most pro- 



nounced at the lower energies and fairly flat over a wide 
reconstructed energy range for the higher true energies. 
This result agrees with the recent analyses by Martini 
et al.[SH] and by Nieves et al. [57] which are based on a 
microscopical model of the matrix element; we take this 
as a further indication that our model contains the most 
relevant physics for this reconstruction. 

While we have performed such analyses only for CC 
processes we note that the formalism can also be applied 
to NC reactions. For the latter a good description of the 
final state particles is absolutely essential for the energy 
and Q 2 determination. 

Since events involving pion production and A excita- 
tions are closely entangled with true QE scattering and 
2p-2h processes, any comparison with data must describe 
at least these reaction mechanisms equally well. In all the 
calculations of the 2p-2h process reported in this paper 
we have always assumed that the outgoing two particles 
from the first, initial interaction are two nucleons. In 
principle there could, however, also be events such as 
v + N% + N2 — > V + N3 + A, or even events with two A's 
in the outgoing state. Such processes would contribute 
to pion production and enhance the pion yield. A solu- 
tion to the puzzle of high pion production cross sections 
mentioned in the introduction might thus come from the 
inclusion of 2p-2h events with one or two A in the final 
state. We speculate that this might solve the MiniBooNE 
"pion problem" shortly mentioned in the Introduction. 
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